EPJ manuscript No. 

(will be inserted by the editor) 



Entropy Production and Particle Yields in Heavy Ion Collisions 
at LHC 



Bela Lukacs and Andras Ster 

MTA KFKI WFI, H-1525 Budapest 114, PO Box 49, Hungary 



(N 

O 

(N 

> 
O 

:z;: 



a; 

(D. 



Received: date / Revised version: date 

Abstract. A sclieme is proposed for calculating the entropy production in highly transparent colliding 
systems. The formalism is continuum physical and fully conform with thermodynamics; and the calculation 
is explicit in entropy production. We have analyzed high energy heavy ion collisions from viewpoint of 
thermodynamics. Entropy densities are calculated for reactions up 5.52 TeV/nucleon-pair energies. Final 
state particle productions are predicted using particle generating quark and hadron models. We find that 
at LHC energy the entropy density is 65.7. From this value we could predict particle yields. The results 
show that there is not yet asymptotic freedom. 

PACS. 05.70.-a Thermodynamics - 12.38.Mh Quark-gluon plasma - 25.75.-q Relativistic heavy-ion col- 
lisons - 25.75.Nq Quark deconfinement, quark-gluon plasma production, and phase transitions 
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1 Introduction 

A group of physicists in Budapest recognised in 1987 that 
calculations of heavy ion collision events need i) thermody- 
namic approach because of the lack of any a priori knowl- 
edge about rehadronisation processes; and ii) the develop- 
ment of a self-consistent anisotropic continuum dynamics 
proper for numeric calculations, because of transparency 

[iiiiiiisjiTiin]. 

Such an approach can be used only once before the 
start of a new accelerator, of course, and then maybe once 
more when the first well-established yield ratios have been 
known, to nail down the actual rehadronisation model at 
that energy. 

In this paper we want to give a prediction for the 
hadronic yields at LHC in Pb+Pb collisions with maxi- 
mal center of mass energy of 5.52 TeV/nucleon-pair. The 
apparatus will work at full energy and intensity in 2014; 
and chances are substantial that even our predictions may 
turn out to be wrong. However, in this case that would 
help, as well. 

The hadronic yields depend on many factors. We know 
some of them fairly well. Some we believe to know fairly 
well. And we are sure that we do not know some well. For 
example, we more or less think that we can predict the 
specific entropy (S/N) at maximal temperature well. We 
believe that the state of Quark-Gluon Plasma (QGP) is 
described correctly. However, we know that we cannot se- 
lect the correct rehadronisation process without empirical 
experience. 
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In any case, we do not yet have this experience at 5.52 
TeV/nucleon-pair. So the strategy will be as follows. We 
use: 

1) A dynamic calculation up to maximal temperature, 
so obtaining the specific entropy. 

2) Then we switch to QGP, exploiting thermodynam- 
ics. 

3) The blob starts to expand and cool, so rehadroni- 
sation comes, for which we use 10 different models, -I- one 
calculation without QGP phase. 

4) The rehadronisation models give various yields. For 
SPS and RHIC energies (158 & 200 AGeV fixed target, 130 
& 200 AGeV storage ring) they can be compared to mea- 
surements, whence we extrapolate for the best rehadro- 
nisation model at LHC energy (at least within the space 
spanned by the present models). 

The rehadronisation models used in this paper are 
listed and briefly explained in Appendix [Xj 

Thus we will give our most preferred prediction. Also 
we will give a few other ones for cases when some factor 
went against our best efforts. E.g. the dynamic calculation 
may be proven invalid (e.g. because of deviation from the 
extrapolation of the differential cross sections). If so the 
prediction is given for another alternative specific entropy, 
as well; and so on. 

This way, if our present most preferred prediction proves 
approximately correct in some months, then we can con- 
clude that nothing unexpected happened. If, in contrast, 
one of the alternative predictions is better, then we at 
least can guess what unexpected happened. Of course, it 
is possible that neither our preferred predictions nor any 
of the alternative ones will be good; then we will have 
been proven brave but unsuccessful. 
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In Section [5] we give a bird's eye view of thermody- 
namic rehadronisation models, while Section [3] calculates 
the specific entropy S/N produced in the collision at SPS, 
RHIC and LHC energies, with flows highly anisotropic be- 
cause of the substantial transparency. Section 2] gives the 
respective numeric values for S/N. 

In Section [5] we select a finite set of rehadronisation 
models from among the infinite possibilities. Which of 
them 8 will be based on the probabilities of quark encoun- 
ters, with or without, however, the inclusions of compress- 
ibility, influence of masses of resulting hadrons and/or 
gluon fragmentation. One more employs sequential fission 
as rehadronisation process, and the last one, similar to 
Model 7, will be defined in due course. In addition, for 
comparison, we will use a model without QGP, as well. 

Of course, we cannot know a priori, which model is the 
best approximation, moreover, this may depend on the 
accelerator energy, as well. So we calculate the hadronic 
yields for the previous SPS & RHIC energies, and then we 
can get an educated guess for LHC. 

This way we select our best candidate, together with 
2 less probable ones (and we repeat the calculations for 
the absence of QGP, for comparison). We definitely do 
have a unique " best candidate" ; the success of one of the 
alternates would mean that some unexpected happened 
between 1 and 5.52 TeV. 

Section[6]discusses some simplifications we deliberately 
have done, and tries to give corrections. Section[7]gives the 
conclusions. 

Some of the heavy ion physicists in Budapst from time 
to time give predictions for future accelerators or thermo- 
dynamic comments to preliminary results back to 1987; 
references will come in due course. The tasks are taken 
with gaps of many years; this is the first prediction from 
us for LHC, assuming that the thermodynamic approach 
works at least moderately well in this regime, too. 

2 On the history of thermodynamic 
predictions of rehadronisation 

Since free quarks cannot be produced, heavy ion physics 
always was regarded as the only experimental way to pro- 
duce Quark-Gluon Plasma. Back to 1979 the predictions 
were that for this either several normal nuclear densities or 
high temperatures (cca. > 160 MeV) are needed. Because 
heavy ion collisions produce high density via collisions, 
heating up is automatic. Therefore, generally predictions 
tended to be optimistic: roughly always: "just the next 
accelerator will produce QGP" . 

Of course, the Budapest group was not immune from 
this optimism. In the middle of the '80's, AGS being con- 
structed at BNL we guessed that 9 GeV/A beam energy 
(on fixed target) would be enough to get pure QGP \1. 
We were not alone with this expectation. True, in 1987 the 
thermodynamic investigations drew attention to the fact 
that the slow decay of longitudinal momentum may cause 
a low-density start of the phase transition |2j which, how- 
ever, cannot end for a while. AGS started in 1986, and af- 
ter preliminary results we believed in QGP at 14.5 GeV/A 



Si-|-Pb [3]; we believed the high A'+/7r+ ratio a signal that 
the K'^ came from QGP. The subsequent data confirmed 
the high ratio [18 , but a methodical comparative analy- 
sis confronting QGP & Hadronic Matter (HM) scenarios, 
including dynamic calculations of both collision and new 
phase formation would have been necessary but then par- 
tially inavailable. For example, we performed a calculation 
taking finite nucleation time in the phase transition into 
account [1], but we did not know, of course, the nucle- 
ation time; and the highly anisotropic local state made hy- 
drodynamics unfounded. So our next analysis could have 
gone only to the breakup we only concluded about a 
quite normal hadronic breakup, with unknown prehistory. 
It became then obvious that dynamic calculations of the 
collisions are inevitable: for this we elaborated the equa- 
tions [51[7], and performed some hydrodynamic calcula- 
tions for the anisotropic matter formed in the collision |5] . 
Since the result was a very extensive transparency pre- 
venting too much density increase, the QGP formation 
seemed more and more improbable at 9 GeV/A (fixed 
target). 

Here it is necessary to stop for a moment. Hadrons 
cannot carry information about QGP; leptons in principle 
can, but for them the conclusions are not easy. We detect 
the hadrons from the Hadronic Matter, after rehadroni- 
sation. Now, the turbulent and unfamiliar state of mat- 
ter at maximal compression may behave as Quark Phase 
from one aspect, and as Hadronic Phase from another. 
Also, very probably, the phase transition is of first order, 
and then, even in equilibrium, mixed phase is possible, in 
which case we may not recognise what happened really. It 
is no surprise that contradictory conclusions were drawn 
from AGS experiments. 

From 1991 attention turned to the SPS under construc- 
tion, at energies more than 1 order of magnitude higher 
then AGS (and to the future RHIC). It was not clear if the 
higher energy would be enough to get QGP; therefore, we 
performed a multimodel prediction, with hadronic matter 
phase as well as QGP with sequential fission hadronisa- 
tion and with various nonequilibrium phase transitions, as 
well. (This model family will be used in the present paper 
too.) For the scheme see [9l[Tgi[TT l [T l [13] . 

When the SPS started to give data, worldwide expec- 
tations were equivocal, but somewhat in favour of QGP. 
However, after 1997 it was generally believed that no well 
developed QGP was present (see e.g. [T^IT^I^ V even if 
some doubts lingered. In 1999 the model ALCOR (Alge- 
braic Coalescence Rehadronisation) was elaborated [14) . 
where quarks are deliberated but do not form a genuine 
phase. While ALCOR itself is not explicitly part of our 
multimodel scenario (being the scenario formed back in 
1991), its philosophy is similar to Sequential Fission, which 
is our Model 10. ALCOR was able to reproduce some SPS 
yield ratios, and we think now that it would have been 
successful even earlier to explain the high if + yields &c., 
unexplained in pure Hadronic Matter scenarios. 

Then came the RHIC in 2000-1, and practically ev- 
erybody expected QGP. However, [T3] was a warning that 
from yield ratios the signature may be unclear. And here 
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observe that we simply cannot make any thermodynami- 
cally "well-founded statement about the transient presence 
of a QGP phase at maximal temperature. For that we 
should prepare and observe systems with different choices 
of independent extensives [H], now particle components, 
so deciding if the prepared systems are similar to each 
other, but if one cannot do this, then the programme 
cannot be realised. The liberated quarks guessed in 1999 
in [14| may quite have been present, but we cannot observe 
it. However, we can observe the hadronic yields. 

From that time (2003) there is a study comparing var- 
ious scenarios when combination/recombination happens 
before detection, so we cannot observe the combination 
processes in themselves. Two of the cases were primor- 
dial hadronisation and heavy ion rehadronisation, for the 
second Models 1-10 were discussed upto S/N=50 [15] . 

We are before the start of data about the yields at 
7000-1-7000 AGeV storage ring energies for p-|-p and at the 
corresponding 2760+2760 AGeV for Pb-^Pb colhsions. No 
doubt, the collision energies will be high enough for QGP. 
However, extreme transparency means that high densi- 
ties are not expected; high temperatures may be present, 
but to decide this the collision must be followed up in 
some anisotropic continuum scenario. As for signatures, 
no doubt, at the same temperature QGP may and would 
result in yields very different from pure Hadronic Matter; 
but the yields are expected to depend strongly not only 
on the transient presence of QGP but also on the actually 
working rehadronisation process; we do not have a priori 
information about the process, and practically the only 
way to decide the dominant way of the rehadronisation is 
to deduce it from the yields. However, we think that now 
we are in a position to make definite predictions. 

Namely, it is almost consensus that QGP was present 
already in the RHIC experiments. So some of our re- 
hadronisation models with QGP, or at least some com- 
binations of them, must give tolerable yields at 65-1-65 & 
100+100 AGeV, and then we can guess the best model for 
2760+2760 AGeV. This will be our method in this paper. 



3 Formulae of entropy production 
3.1 On the anisotropic state 

For simplicity throughout the whole paper we consider 
a matter composed of a single, conserved particle com- 
ponent. This restriction is not necessary at all, could be 
lifted, but now strongly simplifies the equations. This con- 
served component is, of course, the baryons. 

This only component has a flux vector, n*. Except for 
very exotic, probably unphysical situations this flux vector 
is timelike (a spacelike flux would mean acausal motion 
for the particles, and a lightlike one would be a limiting 
situation in the best case). It is algebraically clear that a 
timelike flux vector field always can be decomposed as 

= nw*; UrU^ = 1 (1) 



where is the velocity field of the continuum and n is the 
comoving particle number density. (According to the Ein- 
stein convention of General Relativity, if an index appears 
pairwise, above and below, then summation is automati- 
cally meant.) In local equilibrium the velocity field is 
not only the average of the individual particle velocities, 
but also a significant part of the particles have velocities 
near to m*. However, in the present case of very energetic 
collisions the two original nuclei do not stop each other but 
they interpenetrate. Then the resulted momentum distri- 
bution is not even similar to a "thermal" one, but the 
momenta bins about m * are almost depopulated. This 
is a nontrivial, challenging situation for Thermodynamics. 

Almost any thermodynamic description is possible in 
2 conventions. Either energy E is the thermodynamic po- 
tential, the function of the (independent) extensives, of 
which entropy S is one, or S is the potential, and E is one 
of the extensives. The two conventions are almost equiva- 
lent, and again the exotic details have nothing to do with 
anisotropy. So just for definiteness' sake, let the potential 
be S. 

Then S must be a unique function of the set of indepen- 
dent extensives, moreover a homogeneous linear function, 

S^YbX''^S,rX'' (2) 

(beware the Einstein convention). Without (2) there would 
be no Thermodynamics at all. 

Now, which are the independent extensives? We can- 
not know the answer a priori |21j . Indeed, there are cases 
when even the number of independent particle compo- 
nents is not obvious [16 . While that is 1 in this Section, 
the correct way is always first to take a set of extensives 
which seems to be minimal; and if that is not sufficient, 
we introduce one more and so on. 

In Ref. [2] we manufactured the simplest extension 
for anisotropic local states, with an extra extensive Q. In 
the simplest case Q has a dimension momentum*volume, 
so Q/V is momentum (or, alternately, Q is the momen- 
tum*particle number), and if Q=0 then the local state is 
isotropic. So for order of magnitude we expect Q/N to be 
similar to the measure of momentum anisotropy. 

At this point we chose simple nuclear and quark equa- 
tions of state, and we used some Fundamental Laws of 
Thermodynamics and so. The result was: at anisotropy 
the deconfinement phase transition starts at 5 times the 
normal nuclear densities and ends at 11 times the nor- 
mal nuclear densities. In between there is an n-q mixture. 
With increasing anisotropy the phase transition starts at 
lower and lower densities, but ends at higher and higher 
ones. 

Since we do not know how big the anisotropy will be at 
full overlap, generally we cannot avoid detailed dynamic 
calculations, but the approach does not need the calcula- 
tion of distributions. Anyway, as transparency increases, 
anisotropy increases as well, so nobody should be much 
surprised if any possible experiment would end with n-q 
mixtures. 

Ref. [5] used the simplest nuclear equations of state, 
and a perturbative QCD one for the quark plasma. In |17] 
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we made the anisotropic equations of state for Skyrme 
potential; not as if we believed it but as an explicitly 
momentum-dependent one. If there is a favourite potential 
+ nice arguments, we can repeat the construction in |T7] 
for that one. Now let us focus on entropy production. 

3.2 Dynamics 

We have the more complicated task: an anisotropic, strongly 
interacting, transparent system, starting from maximally 
anisotropic but cold state, and the question is (more or 
less): what will be/ was present at maximal density? (The 
observation is at asymptotic outgoing state, so well after 
maximal density). 

But we have 2 undisputable equations and 1 undis- 
putable unequality: 

= (3) 
n!, = (4) 

s-r > (5) 

The derivations of the dynamic+thermodynamic equa- 
tions are given in [B]; here we give only Schlagworter. 
Everything has an energy- momentum tensor T (see GR); 
the system is compact and separated, so Eq. ^ is true. 
Eq. ^ is true for conserved particles; if they are not con- 
served, write a rhs. And Eq. ^ is the Second Law of 
Thermodynamics. Of course, Eq. ([5|) is an unequality as 
an identity: not just as an accident. These choices of equa- 
tions of state (with some simplifications) will be enough 
as it will be seen. 

We have some anisotropy. It has a direction (at the be- 
ginning surely the beam direction), and a measure. So it 
is a vector. We dub it (anyways, letters tSzu are neigh- 
bours). Surely t* is spacelike. So the local state has at least 
two preferred vectors: timelike unit u*, and spacelike non- 
unit t'. The most general T*'^ with two preferred vector 
fields is: 

+ {bH'' + fb'') + c''' (6) 

where 

drU^ — drt^ — brU'^ — brf — CipU'' = Cirf — (7) 

Vector d^ is well known as heat current; in many cases 
local equilibrium approach neglects it as, e. g., the perfect 
fluid approximation. If so, we may neglect 5', as well. We 
may assume that (a 2*2 tensor) is "as isotropic as 
possible". Then 

T''' :^eu'u^ + f3{uH''+eu'') + k{g'^+u'u^-tH^/f} (8) 

where we wrote a — >■ e due to the usual definition 

T'^u.uk = e (9) 

OK, this seems to be a toy model, and maybe it is too 
simplified. But now the energy-momentum tensor is inher- 
ently anisotropic, and still we shall get in due course all the 



evolution equations from Eqs. (|3l4p and from Uneq. ([S]); 
excepting some simple " material equations" . 

As for the scalar coefficients , equations of state are 
needed; /3 is in a mirror-symmetric state, see [7] for the 
statistical physical approximation; here we remain at /3=0 
for simplicity. The entropy current vector cannot be any- 
thing else than 

^ su' + zf (10) 

because of the approximation (7). As for the coefficient z 
we shall get an equation immediately. 

Then Eq. ([3]) is an evolution equation for the extensive 
density e, plus the equation for the acceleration (so for u*). 
Eq. (lU is the equation for n. We still need an evolution 
equation for the anisotropy , but we have Uneq. ([5]). 

Here a short interlude is needed. Obviously the extent 
of anisotropy is characterized by the absolute value of the 
anisotropy vector; this quantity will be called t, which is 
now not the time. (Sorry for the convention; it goes back 
to 1990.) We may use the Ansatz that the direction of is 
the beam direction. However, surely, t is not an extensive; 
it looks like rather as a specific extensive. However, a good 
candidate seems to be e.g. 

Q = Vnt (11) 

and then the extra extensive density is 

q = nt (12) 

and the proper thermodynamic potential density s is 

s = s(e, n, q) (13) 

in the simplest case. However, just for now for simplicity 
let us use the noncanonical quantity s 

s = s(e, n, t) = s(e, n, q = nt) (14) 

For a straightforward but tiresome step see Appendix |B] 
The results is a few equations and unequalities as follows: 

k — p{e, n, t) -f ujTs,t +Su^^^ 
p = T{s - ns,n -es,e ) 
z^/3/T 

U!>0 

S>0 

l/T = s,e (15) 

where 

Dt= X + l^{Tr + TUr:sU')f + OffUr-.s + WM^,. (16) 

and 

s,t A > 

e = ik/f~j)/s,t (17) 

Eq. (jB.ip shows that A contains the characteristic time 
for the decay of the longitudinal (=beam direction) mo- 
mentum, so a weighted cross section. As we told, just now 
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/3 = 0. As well, we can take the approximation uj = S = 
for the present. And then A is the only new quantity 
needed. It governs "thermalisation" . 

Then we have got the full system of evolution equations 

as: 

Dji + nur-r = 

De + {e+p + qju'.^ = 

Du + {D{p + q) + {p + q)„}/{e+p + q)^Q 

Dq + {q + q/v)u^^- X = Q (18) 

where x is the beam direction, and for the connection of 
q and its canonical conjugate v see [T7]; the connection, 
of course, depends on the equations of state. The system 
(Eq. (US])) contains the evolution equations for all inde- 
pendent extensive densities (now for simplicity's sake e is 
one of the extensive densities and s is the potential den- 
sity; remember Eq. (IB.1[) ') plus for the velocity; so we can 
calculate the local state plus the energy-momentum tensor 
the next moment. 

We are now ready with the evolution equations; still 
the equation of state must be specified to establish the 
extensive-intensive relations, and also A, the coefficient of 
the decay of longitudinal momentum must be calculated. 

The quantity v is the new entropic intensive canoni- 
cally conjugate to Q=Vq. Its actual functional form de- 
pends on the equation of state. In [S] and [7] we used 
a simple enough model when the matter is dominantly 2 
cold interpenetrating particle currents, still there is a tem- 
perature T with a simple parabolic compression potential 
and a gluon radiation field. |17| is nicer, with a Skyrme 
potential in the hadronic phase; but now we are giving 
only the formulae for the quark phase, where 

p = nT+ (7rV30)r'* + {K/9)nox^{v arctg{v) - 1 
+x{2/y~y)) 

e = mny + {3/2)nT + {tt^/10)T'^ + (K/18)n{x - if 
+ {K/9)niy - l){x/y^iiy + l){x/y - 1) 

q = mxnQln{v + y) + (K /9)nQ{{x/2)ln{v + y) 

-x^arctg{v) + vx^{{ix/2y) -Ijy^ - xv^/y^)) 

y^{l + v^Yl^ 
X = n/no 
no = 0.16/m~3 

m = 938MeV (19) 

These equations are somewhat implicit, but hence v can be 
calculated as function of e, n and g, and then the evolution 
Eqs. are complete indeed. 

As for A, Eq. shows that it is the decay rate of 
the density of extra momentum in beam direction nt. In 
ultrarelativistic situations, as now, the momentum distri- 
bution in beam direction has two peaks with widths ~r 
and separation '^t'^Ebeam- So in first approximation the 
distribution is sharp. 

Now consider a collision. Most collisions happen be- 
tween particles moving oppositely with almost ±Pheam- 



So it is enough to evaluate such a collision. The details go 
into Appendix [C] because of unavoidable inconsistencies 
in the notation. 

We would need differential cross sections, but that is 
not easy for inelastic collisions at such high energies. We 
chose the compromise to use the total (elastic -|-inelastic) 
cross section with the shape of the elastic one in the mo- 
mentum transfer variable. For the latter we used the curve 
of Islam and al. [H] . 

3.3 Some results and outlook 

We solved the dynamic equations in the range of E/A ^ 1 
GeV long time ago [8J . But the equations can be solved in 
the 100 GeV - 10 TeV range as easily as below. Ref. [7] (at 
1 GeV/ A) gave quite rational results (for 1 GeV/ A). The 
particle density is cca. 2no at the center, cca. no a-t the pe- 
riphery and evolves in between. At the beginning T is cca. 
150 MeV in the center, dropping towards the periphery; 
later slightly higher well inside. The anisotropy is highest 
somewhere between the center and periphery, decreasing 
slowly in time. Indeed, anisotropy exists when both left- 
moving and right- moving particles are present in the same 
volume element; because of transparency some anisotropy 
remains even until the end of the overlap (and afterwards 
the calculation must stop, with a stage of breakup). 

At 200 GeV/nucleon-pair the transparency will be higher. 
So while the kinetic energy would be much enough for 
deconfinement phase transition, it would be transformed 
very ineffectively into temperature, and even more inef- 
fectively into compression. So without full dynamic calcu- 
lations there is no unequivocal prediction for or against 
deconfinement phase transition. 

But the full dynamic system of equations is given above 
and can be used. And such a scheme is explicit for/in the 
generation of the entropy. 



4 Numeric results 

Our simplifications for the entropy production of the col- 
lision are as follows: 

a) Slab geometry (it does not seem too important). 

b) Homogeneity (more or less true due to interpenetra- 
tion). 

c) The equations of state given in Section [3] with K=0. 

d) No collision after crossing over (this is more or less 
true). 

e) Breakup just afterwards. 

With these simplifications and with Islam's differential 
cross sections [22] we have numerically solved Eqs. (|T8]) 
and we have calculated the S/N production for the accel- 
erator energies used so far, as well as for LHC. The results 
are shown in Table [T] and Fig. [T] 

We have calculated S/N of central Pb-I-Pb collisions 
for 6 CM energies 17.3 GeV, 19.4 GeV (158 & 200 GeV 
fixed target), 130 GeV, 200 GeV, 7000 GeV & 14000 GeV, 
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Fig. 1. Calculated specific entropies at different center of mass 
energies. The solid line is drawn for guiding eyes. 



ECM [GeV] 


S/N 


17.3 


51.4 


19.4 


52.2 


130 


57.5 


200 


58.6 


2760 


64.2 


5520 


65.7 



Table 1. Specific entropies at various energies of center of 
mass. 



respectively. While the last 2 energies seem to be possible 
only for p+p collisions, the calculations are meaningful, 
as it will be seen. The results are shown by Fig. 1. 

The analysis of the curve shows that a power fit is not 
too bad but with a hardly interpretable power law 0.03, 
while a logarithmic fit S/N w 45.31 + 2.381n(£;cA/ [GeV]) 
is even better. The logarithmic behaviour is not surprising 
in an ultrarelativistic regime without energy scale. 

We cannot determine the error of the calculation, but 
we can give mean deviations from the logarithmic curve. 
However, doing this one can see that the main deviation 
would come from the two SPS data points; obviously be- 
low 20 GeV we are not in the true ultrarelativistic regime. 
Removing these points the fit is much better, resulting in 
S/N « 47.019 + 2.168 \n{EcM [GeV]), with a error 
^(S/N) = 0.14. The error is quite low, and we will use it 
from 200 GeV onward. So, for the relevant energies the 
specific entropies are given in Table 1. The dynamic cal- 
culations give S/N=66.1 for 7000 GeV and 67.8 for 14000 
GeV. 

The entropies are in the expected range. Indeed, Letes- 
sier, Rafelski & Tounsi read out almost these values from 
yields and energy distributions cca. 15 years ago. At Ej^eam 
/ A = 200 GeV, which is the second energy of the Table, 
Vs = 19.4 GeV. In 1993 the value was 50 ±4 [20], in 1995 
42 ( S/N ( 48 [23] and in 1996 35 ( S/N ( 60 [TO], any of 
them an excellent first guess. 



5 Application of rehadronization schemes 

The most direct experimental facts about QCD are the 
particle yields at the detectors. However, these particles 
are hadrons, coming from a hadronic phase after rehadro- 
nisation. Only their quark constituents keep the memory 
of the QGP phase. Therefore, the actual rehadronisation 
process seriously infiuences the yields, but we do not know 
a priori the correct rehadronisation model. 

Here we use a variety of rehadronisation models simul- 
taneously: each one would lead to a set of yields, and at 
the end one of them or at least a combination of them 
might be satisfactorily near to measured yields at lower 
energies. Then we may take this model or the combina- 
tion or an extrapolation. This decision will be done in the 
next Section; now let us see the models used here. Calcu- 
lations in models 1-10 were made a few times; of course, in 
contexts of much lower energies [^ [TUlllll[T^[T5] . We con- 
sider here 10 rehadronisation models (-1-1 without QGP) 
from amongst the virtual infinity of them. 

Models 1-8 assume a QG plasma at maximal com- 
pression; in the QG plasma there is complete thermal, 
chemical and mechanical equilibrium (processes are rapid 
enough). Then an expansion starts and somewhere the 
matter reaches the phase boundary. As a matter of sim- 
plification the phase boundary is fixed on the p-T plane, 
and for the actual choice see jMj. Then rehadronisation 
starts, but the growing HM droplets remain in equilibrium 
with the QG environment (no supercooling &c.). The pro- 
cess goes until the last remnants of QGP vanish: then we 
detect the hadrons. 

Even this scenario means an infinity of models, of which 
here we take 8, as follows. We have 3 properties of the 
models, and 2^ = 8. 

1) Gluon fragmentation may occur during rehadroni- 
sation: we may neglect it (0) or may calculate it as 15% 
ss and 85% qq. For the suggestion of this particular ratio 
see [H]; for most probable strength see PBI12II24| . 

2) The probability of hadronic combinations may come 
from the quark numbers by simple combinatorics (0) or 
they may contain an extra factor ^ e^^l^ (1), where F is 
the free energy. This latter idea is obtained from the num- 
bers of final state microstates ^ [22 , and for technical 
reasons here we will approximate F by the hadronic rest 
mass m. 

3) The hadronic phase may have a compressibility of 
2AQMeV/fw? (1) or not (0). 

The 3 numbers are regarded digital, and then add 1: 
so they are indeed Models 1-8. Model 9 is the absence 
of QGP phase at all, for simplicity's sake with complete 
thermal, chemical and mechanical equilibrium throughout 
the whole collision, and Model 10 is a sequential fission 
one [TT]. A comparison of the results of Models 1-10 was 
given in [TS] up to SPS energies. We will add Model 11 in 
due course. All of these hadronisation models end in sta- 
bles. Models ending in resonances are not considered here, 
because very various such models would be possible. We 
will later discuss the consequances of this approximation. 
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Fig. 2. S jn ratios of model predictions above S/N=5. The 
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Fig. 3. H/n ratios of model predictions between S/N=0 and 
10. The solid lines are drawn for guiding eyes. 



12 hadrons are included into the scheme, with the as- 
sumed QG precursors, as indicated: 

N = qqq 

Y{= AhJ:) = qqs 
S — qss 
f2 = sss 

+antibaryons + 
TT^qq 
K = qs 
K = sq 

T]^ss (20) 

for a while without distinguishing u and d. All of these are 
"stables" , i.e., they do not decay in strong interaction. Ob- 
viously in the initial state there is a slight d excess, which 
may cause a bias (see next Section), but not too much, 
and the last mesonic representation is highly debatable 
but rj will not in fact be used explicitly. 

Then we get unique predictions for the 12 yields if the 
model and the specific entropy are fixed. To demonstrate 
the model and specific entropy dependences, consider the 
ratio E!/n, for the Models 1-10 between S/N = 5 and 95 
and separately between and 10 on Fig. [5] and Fig. |31 
respectively. 

As it can be seen, for Models 1-8 there are 8 different 
predictions at low S/N's, but model pairs practically co- 
incide for S/N>20. So henceforth, it will be superfluous 
at RHIC and LHC energies to refer separately to Models 
2,4,6 & 8. 

Of course, very probably none of the rehadronisation 
models will be satisfactory because i) the rehadronisa- 
tion processes do not clearly separate; and ii) some of the 
present models have fixed parameters which are not nec- 
essarily the true ones. Problem i) is seldom formulated ex- 
plicitely, although it is serious. As an example, it is quite 
possible that in rehadronisation the result is not p = uud 



but Z\+ = uud with parallel spins. Then Z\+ decays in the 
hadronic phase giving extra mesons, mainly pious. None 
of the the 11 rehadronisation models include resonances, 
for obvious difficulties. One of them is that very probably 
some uud triads would correct the spin positions still in 
quark phase, while some Z\+s are formed and then decay 
in the hadronic phase; and we do not know, how much. 
(There is some indication for such decays and extra tt's 
at 200 AGcV RHIC as it will be seen in Section [5]) We 
think only tt numbers will be sensitive on this. As for to 
overcome the second problem we can try with linear com- 
binations of models with weights fitted to the observed 
yields. However, at 200 GeV the available yield ratios are 
limited. So we give here the method only for 3 models se- 
lected in the same time. Let the models be denoted by f, 
g and h, respectively; the index i denotes the actual ratio. 
Then the actual yields are given by the combinations 

i?, = a{E)f, + b{E)g, + (1 - a ~ h)h, + a, (21) 

where the cr^'s are the errors of the combined model, to be 
minimized in average. This means that 

{N - 2)a^ = Sf^^af = min. (22) 

Hence 

a^iIi(/» - h,)' + bS^^.if, - h,){g, - h,) 

= St^{R,-hi){},-hi) 

aStiih - h^){g^ - h,) + hI:tl{9^ " h^f 

= S^^iiR, ~ h,){g, - hi) (23) 

Hence we get the best weights, and substitute them 
back Eq. ([^ gives the average . 

Now we are going to compare the measured particle 
ratios at 158 & 200 GeV fixed target (SPS) and 130 & 200 
GeV storage ring (RHIC) experiments to the predictions 
of the rehadronisation models 1-8, 10 plus Model 9 without 
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Model 


158 GeV 
SPS 


200 GeV 
SPS 


130 GeV 
RHIC 


200 GeV 
RHIC 


1 & 2 


0.394 


0.759 


0.012 


0.0503 


3 & 4 


1.326 


2.819 


0.009 


0.0306 


5 & 6 


0.120 


0.141 


0.005 


0.0330 


7 & 8 


0.165 


0.283 


0.005 


0.0308 


9 


0.020 


0.028 


0.072 


0.0950 


10 


0.084 


0.176 


0.054 


0.1288 



Table 2. The average predictive squared errors of the indi- 
vidual models. 



QGP to find the best one or the best combination. We 
took the ratios from [TlllinilllliaiSZlllSlllSllSniEIlElISS]; 

when we found more than one ratio but still near to each 
other, at a few a, we averaged with least squares. (When 
they were too different we discarded one of them.) We 
tried to use low impact factor collisions, as far as possible. 
Later Table [7] will give the ratios we used in the analysis. 
The deviation a is the quantity in Eq. ^F^: the average 
"error" of the rehadronisation model. It is dimensionless. 
Table [2] gives the average squared errors for individual 
models. However, note that if a single model is evaluated 
the divisor changes from N-2 to N-1. 

158 GeV fixed target means 17.3 GeV in storage ring 
situation, at least for symmetric situations. For 200 GeV 
the equivalent energy is 19.4 GeV. As shown by Fig. 2 at 
the entropies relevant now the first 8 models give only 4 
different results up to 3 digits. 

The picture is clear enough: at the two lowermost ener- 
gies the best predictions are given by Model 9, so a model 
with no QGP at all. Even Model 10, a sequential fission, is 
much worse. The higher error at the higher energy prob- 
ably comes simply from the fact that we used less ratios 
there; still Model 9 wins clearly. So it seems that at those 
energies no well developed QGP was formed. This seems 
to be conform with the opinion of Ref. [H] and [31] from 
1996-7 even if in the previous years still the hopes were 
high for 19.4 GeV, and even if Ref. [34] in 1999 still was 
optimistic. Just oppositely, for 130 & 200 GeVs Ref. [55] 
expects QGP in 2006, and indeed at this energies QGP 
models win. At 130 GeV the best models are 5 through 
8; at 200 GeV even 3 & 4 are as good, but at the highest 
energy we could use only 5 ratios. So we are giving greater 
weight to the 130 GeV errors; but clearly we have to try 
with the linear combinations, as well. 

There is no doubt at the two lower energies. The three- 
component combinations of smaller af always contain Mod- 
els 9 or 10; if both then the weight of 10 is either smaller or 
even negative. With a dual combination of 9 & 10 the error 
minimum is at combinations with weight slightly above 
1 for Model 9 and slightly negative for Model 10. Since 
weights either ) 1 or ( are aphysical, this means that in 
the whole 9-dimensional space of the linear combinations 
of Models 1-10 the physical optimum is at (or very near 
to) the pure Model 9. 

For RHIC at 130 GeV the situation is more compli- 
cated. Models 5-6 and 7-8 give very similar small errors. 
Trial combinations with both 5 and 7 give various min- 



ima but always at aphysical weights. Dual combination of 
5 and 7 is optimal at positive weight of 7 and negative 
one for 5. Finally an analysis of 5-7 combinations gives 
monotonous decrease in the physical domain from pure 5 
to pure 7, but the whole difference is small: for Model 5 
crf=0.0061, for Model 7 0.0058. 

At 200 GeV the result is similar, although the overall 
error is higher (less ratios were used) and Model 3 is very 
sligthly better then Model 7. The squared average is gen- 
erally high (which is the reason that we remain at Model 
7 as best). We refer to this point in due course. 

All the models 5-8 contain gluon fragmentation. So at 
least from the particle yields, it seems as if both a well- 
developed QGP and substantial gluon fragmentation ap- 
peared somewhere in the range between 20 and 130 GeV. 
It would be interesting to see yields between 20 and 130 
GeV; but for now we go directly to 5520 GeV. We cannot 
expect Models 9 & 10 to overcome again at high entropy; 
our best educated guess is Model 7(-8). However, we give 
predictions of Model 5, as well, having been almost as 
good at 130 GeV. Because of the substantial presence of 
gluon fragmentation we give a Model 11, as well, which is 
Model 7 but with exaggerated gluon fragmentation, tuned 
up to the energetically possible limit. (In our language it 
is 1.36 instead of 1 for the gluon fragmentation param- 
eter [I(71ll2) .) Also, Model 11 may mimic somewhat re- 
hadronisation into resonances (both processes resulting in 
extra light mesons). 

Now, we are in the position to give the 5.52 TeV pre- 
dictions. (This energy is valid only for Pb+Pb, being then 
Z=82 and for average isotopic coposition A=207. How- 
ever this energy is almost constant in a wide range, e.g. 
in Au+Au it is 5.61 TeV, for example.) First we give Ta- 
ble|2]for stables defined in Eq. (j20p . still not distinguishing 
between the light quarks. Then we list 4 simplifications 
inherent in the calculations but easy to estimate for the 
errors caused, and even correct it in some cases. 

Point 1. Slight Simplification: Insensivity on Charge 
and its Consequences 

So far, we have not distinguished between the 2 light 
quarks. Therefore, N(p)=N(n) &c. This is definitely not 
true for the initial state with heavy nuclei, although it 
is still not an absurd approximation. The final results 
can somewhat be corrected according to the u/d ratio 
of the initial condition; but then the ratios should be 
calculated individually for different beams. For Pb-fPb 
N(M)/N((i)=0.87. Here we are going to make the correc- 
tion. 

In the charge insensitive approximation: 

E N{u) = N{d), then 
N{p) = N{n) 

N{U+) = N{U-) = N{U") = N{A)/2 

and similarly for antibaryons, 
N{tt+) = iV(7r") = N{tt°)/2 
N{K°) = N{K-) 

and similarly for the ant ip articles. (24) 
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Particle 


Model 5 


Model 7 


Model 9 


Model 11 


Model 7 
(with S/N=75) 


Mass 
(MeV) 


N 


528.79 


475.95 


389.19 


585.18 


512.60 


939 


iV 


326.50 


301.58 


152.70 


365.80 


350.83 


939 


Y 


610.63 


711.38 


261.25 


658.50 


801.15 


1174 


Y 


412.91 


479.61 


113.76 


461.67 


557.46 


1174 




79.09 


38.49 


59.65 


26.79 


44.11 


1318 




64.00 


30.98 


32.14 


22.33 


35.94 


1318 


n 


19.79 


1.01 


8.03 


0.35 


1.17 


1672 


n 


19.79 


1.01 


5.53 


0.35 


1.17 


1672 


TV 


3755.23 


3686.70 


1775.05 


4104.40 


4212.90 


138 


K 


1285.38 


1388.59 


965.63 


1273.83 


1579.55 


496 


K 


1057.50 


1141.80 


757.30 


1068.06 


1319.53 


496 


V 


162.89 


125.60 


358.21 


96.23 


144.64 


549 



Table 3. The quantities of Point 4 at 5.52 TeV. Y stands for {A, S). The bold column is our favoured prediction. 



(For the Y's see App. C.) 

Now, consider a well-developed QGP. Really, N(d) > 
N(u). However, the asymmetry is smaller then in the ini- 
tial state, because the processes creating quarks & anti- 
quarks are charge symmetric. At the state just before the 
rehadronisation 

N{d) + N{u) = 1242 + N(d) + N{u) 

N{s) = N{s) (25) 

(The second formula is invalid for Model 10, and there 
are no quarks at all for Model 9.) 

All antiquarks have been produced in charge-symmetric 
processes and all quarks above the initial 1242 ones, as 
well. Since the rehadronisation uses the same pool for any 
hadron, the final result is (with R = d/u) cca. as follows. 
Left hand sides are the charge asymmetric yields, right 
hand sides are the quantities in Eq. (^0)) . 

p = N/{R+l), p = N/2 

n = N^R/{R+l), n = N/2 

yl = r * /2, A = Y/2 

S- = Y ^ R^/2{R+lf, Z^ = F/8 

Z'° = y*i?/(i?+l)2, Z^ = F/4 

S+ = Y *l/2{R+lf, r+=F/8 

S° = X *1/{R+1), S0=X/2 

S- = X * R/{R+l), 'W^^X/2 

n- ^ Q, JF = 12 

7r+ =7r*(l/2)/(i?+l), 7r°=7r/2, tt" 

= 7r* {1/2)* R/{R + I) 

K+ = K *1/{R + 1), 'K+ = K/2 

K" = K * R/{R + 1), 'K" = K/2 

rj = 7] (26) 

For colliding nuclei others than Pb the charge-symmetric 
Table [3] scales with Nb, the initial u+d = 3Nb, and R can 
be calculated from the initial u/d and from the rescaled 
Table El For Au-I-Au N(u)/N(d) = 0.88. 



Point 2. A Minor Neglection in the (S,A) Segment 
There is a slight difference between the E and A masses. 
This here will be ignored (as it will be seen, is not really 
observable, anyway). There is an even slighter difference 
amongst the S masses, also ignored. The reason of the 
mass differences is among Ss and between H's that the 
quarks are the final components. So minimal S energy is in 
uus combination. These differences are a few MeV's. The 
mass difference between A and is, however, 77 MeV. 
While this may still be not too serious, the real prob- 
lem is that both have uds composition the relative spin 
positions differ. Therefore, in the {A, S) sector the reso- 
lution depends in somewhat arbitrary assumptions until 
resonances will be included. Our present assumption is 
discussed in Appendix |D] it gives the numbers in Eq. (I26p 
and this problem will also appear in next Section. 
Point 3. The Stables 

Our model includes the "stable hadrons" of the parti- 
cle physicist convention, that is those which would be sta- 
ble were the s quark stable (no weak interaction). These 
"stables" are: 

E+,S\E-,A,S°,S-,Q- 
the respective antihyperons, and the mesons 
TT+y,TT-,K-,Kl,Kl,K+,i^. (27) 

They do not decay on 10~^^s scale. Without distin- 
guishing between w an d (charge insensitive case) there 
were 12 stables as in Eq. (1^0)) : distinguishing them there 
will be 26 ones as in Eq. (j26D . However, really r] is rather 
obscure and it contains both qq and ss. Now, observe that 
17°, 7r° and 77 are not detected. (They are "stables", but 
they do not reach the detectors). The lifetimes are long 
enough to survive breakup but not long enough to di- 
rectly detect them. So it is only of secondary importance 
what is the exact quark content of 7r° & 77. (Really, in 
the nonett two heavy and not too familiar measons, 77 
and 77', are represented by our symbolic 7/.) Here we regard 
them as the lighter composed of a light quark pair and the 
heavier of a strange pair. In addition obviously, Models 1- 
8 contain a very important neglection, mentioned already 
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briefly. Namely, some hadrons probably will be first res- 
onances not stables. But both rehadronisation and reso- 
nance decay are in 10 ^^'^s characteristic time, so some res- 
onances decay in the hadronic matter, giving extra light 
mesons, mostly pions. This process is not included in the 
models, and it would be rather difficult to do it. Here, we 
only declare that the numbers are very probably underes- 
timated, therefore, in all the models. 

Also, there is the well known K'^ problem. For all this, 
see the next Point 4. 

Point 4- The Observables 

For handling 17", 7r° & 77 observe that neglecting lep- 
tons & photons IJ° decays solely to A, tt" decays solely to 
nothing, and 77 gives a 7r° (which then decays into noth- 
ing) or either nothing else or a 7r"'"7r~ pair. Since Z'", tt" 
& r] decay before the detectors, in them we get: 



7] ^ 

7T+ + 0.29 * r] 

TT" ^ TT" 0.29 * 7y (28) 

where — >■ does not mean any temporal sequence but a 
calculation of detectables. For K^^s the detectable com- 
ponents are and Kg, for which 

Kl = K'^^ = {K° + K")/2 (29) 

This way we still have, say, Z'" if it can be reconstructed, 
which is sometimes claimed. 



6 Discussion 

Earlier we have told that at lower energies (158 & 200 
GeV/nucleon fixed target) the model best for reproduc- 
ing the yield ratios amongst the 11 ones included in the 
present survey was Model 9, the only one without QGP 
phase at all. This is a signal that at these energies pure 
quark phase is still absent. However, at 130 & 200 GeV 
per nucleon storage ring energies the two best models are 
5 ("100") and 7 ("110"). Model 7 is marginally better, and 
linear combination does not improve the performance: at 
200 GeV the situation is not clear but Model 7 is better 
then Model 5. Note that individual yields may differ be- 
tween Models 5 & 7; only they are near to each other in 
overall performance. 

Since this fact indicates well-developed QGP already 
at 200 GeV, we do not expect anything else at 5.52 TeV. 
Therefore, our favourite prediction is Model 7; but we do 
not deny a slight chance from Model 5. As for Model 9, 
the one without QGP, we give its yields only for compar- 
ison. Model 11 is Model 7 but with the maximal possible 
gluon fragmentation. Finally, we give also the predictions 
of Model 7 for S/N = 75 (instead of the expected 65.7); 
the reason will come. The final results are in Table HI then 
you can calculate any ratios of the observables taking into 



account Eqs. (^51 (23). As for the error of the prediction, 
in theory the test at 200 GeV suggests ±0.02 in ratios, 
but the (t's in Table [5] show that there are serious biases, 
too. 

Of course most biases/systematic errors are unknown, 
but we can estimate 2 of them. 

First, as told in Point 1 of Section [SJ the simplified 
model overestimates the number of quarks u. Because of 
heavy nuclei in the beam the true u/d ratio is « 0.87 {Pb+ 
Ph). Now, all in the ratios /p, Q~ / S~ and tt^ / K~ the 
ratio N(s)/N(d) appears, and N(d) is in reality greater 
than in the model by cca. 6 % relative. However, note 
that this distortion is less in Model 7 where the gluon 
fragmentation somewhat dilute the d surplus, and even 
less in Model 11. And in Table [H] we will give the charge- 
sensitive yields according to Eq. (pS)) . 

The second type of error may come from the thermo- 
dynamic nature of the models; in them particle number 
conservations are not explicit. But the final particle num- 
bers can check the accuracies. Table |4] summarizes some 
results of Table 121 

We can see that Models 7 & 9 are excellently consis- 
tent with the conservation laws, and the relative errors 
are low even for Models 5 & 11, because the baryon num- 
ber & strangeness conservations should be maintained in 
the presence of lots of particles with negative quantum 
numbers. If we relate the above errors to the averages of 
numbers of particles with -f & - quantum numbers, then 
at the worst cases the relative error is 0.1% for strangeness 
and 0.7% for baryon number. 

Observe the tremendous particle production in the col- 
lisions. Just after breakup there are at least an order of 
magnitude more particles present than in the initial con- 
dition. As for the differences between the total yields of 
the models we see some 1% difference between Models 5 
& 7, hardly more between 7 & 11 and the only serious 
difference is Model 9 without QGP phase. 

Charge conservation will be checked after the charge- 
asymmetric Table ini 

Table [5] gives us the data for estimating the efficiency 
turning the beam energies into " something" . At the initial 
condition there are 414 nucleons, each with 2760 GeV ki- 
netic energy. Just after rehadronisation/breakup some of 
the kinetic energy has turned into rest mass, some into 
internal energy in thermodynamic sense, and some en- 
ergy has gone away in the beams. Calculating with a 160 
MeV temperature at rehadronisation as approximately got 
in Models 1-8 & 11, the numbers are as follows: where 
Etransf docs iucludc the rest masses, but not the ones 
present in the beginning. At the initial condition Etransf 
was 0, and Mtot/414 was 0.94 GeV. 

So we can see that the efficiency of conversion is pre- 
dicted to be (3 - 5) * 10-3 at 5.52 TeV. The reason is of 
course the forward-peaked differential cross section. (As 
told earlier if rehadronisation goes through resonances, 
not included in the models, then the resonances yield ex- 
tra light mesons. This, of course increases the kinetic en- 
ergy transferred to masses. However, the effect is not big. 
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Total 7V° 


Initial 


Model 5 


Model 7 


Model 9 


Model 11 


Particle 




8322.50 


8382.70 


4878.44 


8653.47 


Strangeness 





0.02 


0.00 


1.68 


-0.01 


Baryon 


414 


415.10 


413.65 


413.99 


420.68 



Table 4. Conservation checks for Pb+Pb collisions with 5.52 TeV/nucleon-pair. 



Model # 


£tot/414 [GeV] 


Mtot/414 [GeV] 


Etransf/ 414 [GeV] 


5 


2760 


9.73 


13.62 


7 


2760 


9.80 


13.72 


9 


2760 


5.77 


7.66 


11 


2760 


9.78 


13.86 



Table 5. Conversion of initial kinetic energy Etot into energy transferred to new particles. 



Doubling the number of pions the last column of Table [S] 
goes up with cca. 1.5) 

At last, we are going to investigate the unavoidable 
consequences of the present lack of observations at ener- 
gies between Tevatron and LHC. In Section 2] we evalu- 
ated the expected LHC specific entropy and we got 65.7. 
However, the experimental differential cross sections are 
not yet measured above Tevatron momenta; we have used 
extrapolations [2 2) assuming that nothing new happens 
there. This is probably so; but if we cross the energy range 
of "something new" between Tevatron and LHC energies, 
then in first approximation a "reasonance peak" is ex- 
pected, so increasing the cross section and so the S/N. 
This "something new" might be anything (e.g. crossing 
the supersymmetry energy scale) and we have no preferred 
candidate, at all. It will not be necessary either; we can 
simply perform calculations with our best. Model 7, with 
S/N = 75 (some 10% increase) for the observables accord- 
ing to Point 5. 

Now, we make the charge-sensitive evolution of yields 
for Pb-|-Pb. Accepting the charge-insensitive yields of Ta- 
ble [31 we operate as told earlier for u's and d's in Point 1, 
so finally: 

N{u) = 578 + {l/2)[N{q) - 1242] 
N{d) = 664 + {l/2)[N{q) - 1242] 
N{s) = N{s) 

N{u) = N(d) = {l/2)[N{q) -1242] (30) 

with this we get Table [6l with stables, with our definite 
prediction in column 3, Model 7, S/N=65.7; if any of the 
other columns will have been proven superior, one can at 
least guess hence the reason. For the detected ones one 
can use Eqs. (^51 We have chosen this way, because 
indirect methods may reconstruct at least 17" and 7r° . For 
comparision, the last column gives figures of Model 7 for 
200AGeV at RHIC. At this point we can check charge 
conservation. Interestingly, charge conservation depends 
on the handling of the (A, U) sector. We tried two dif- 
ferent decompositions of K's into A and S which will 
be discussed in Appendix |DJ Exact charge conservation 
would mean Z=164; at the first choise Z « 174 for all the 
models, in the second one Z w 172. So we used the sec- 



ond for yields, but our calculations seem to indicate that 
further supression of Us do not help in charge conserva- 
tion. Note, that our rehadronisation models are thermo- 
dynamic, charge asymmetric states are not of minimal en- 
ergy. Therefore, thermodynamic models have the tendency 
to violate charge conservation. Anyway, the net error 5% 
correspond to error/all charged « 0.2%, in the order of 
the errors for strangeness conservation. 

In Table [7] we give all the abundance ratios which we 
used in Section [5] the data will show the qualities of the 
models in the last column. The experiments and the refer- 
ences are identified, too. We have ignored individual data 
that seem to be direct contradictions to data referred here. 
In case of 158 AGeV SPS we have not used K'^ /tt^ and in 
case of 130 AGeV RHIC we have not used Q/-k~ , /tt~ , 
p/tt~ because we were not convinced that the pion back- 
ground does not adulterate the ratios. (Again, hadroni- 
sation into resonances being neglected) In contrast, such 
ratios were kept for 200 AGeV RHIC because of the short- 
age of ratios at this energy. If one does not like such ratios, 
one can use a derived experimental ratio K~ /p instead of 
the last 2 rows given in the brackets. The predictions of 
the models with least average errors are given, too: Model 
9 for the first 2 energies and Model 7 for RHIC energies. 

Now, let us compare methodically the measured and 
predicted ratios. Generally, they must agree because the 
model errors were estimated in Eq. (|22p and Table [5] from 
the differences. However, it is interesting to see the possi- 
ble individual big differences which would signal shortcom- 
ings of even the best model of the original 1-10. Henceforth 
we study the deviations ri^obs — fi^mod whose errors are 
composed from the individual measurement errors and the 
model error. As for the error of the difference A, its distri- 
bution is the convolution of those of the r's, so the second 
central momentum is the sum of the constituents 136 , i.e. 

a^^al+ai (31) 

but of course we do not know the distribution, only that 
it is nearer to Gaussian than those of the r's. 

Using the good rehadronisation model, the observed 
r's should fluctuate around the model predictions, so the 
set Ai denoting individual ratios should fluctuate around 
0. Table E] contains 26 independent ratios. Which of them 
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Particle 


Model 5 


Model 7 


Model 9 


Model 11 


Model 7 
(with S/N=75) 


Model 7 

(200 GeV RHIC) 


P 


261.53 


235.41 


190.87 


289.62 


253.85 


221.30 


P 


163.25 


150.79 


76.35 


182.90 


175.42 


131.99 


n 


267.26 


240.54 


198.32 


295.56 


258.76 


226.66 


n 


163.25 


150.79 


76.35 


182.90 


175.42 


131.99 


E+ 


74.68 


87.01 


31.42 


80.65 


98.23 


78.45 


E+ 


51.61 


59.95 


14.22 


57.71 


69.68 


52.52 




152.64 


177.82 


65.29 


164.61 


200.27 


160.69 


EO 


103.23 


119.90 


28.44 


115.42 


139.37 


105.04 


E- 


77.99 


90.85 


33.92 


83.99 


102.07 


82.29 


E- 


51.61 


59.95 


14.22 


57.71 


69.68 


52.52 


A 


457.95 


533.51 


195.91 


493.86 


600.84 


482.12 


A 


309.68 


359.71 


85.32 


346.25 


418.10 


315.13 




39.12 


19.04 


29.25 


13.26 


21.84 


16.90 


HO 


32.00 


15.49 


16.07 


11.16 


17.97 


13.60 


S~ 


39.97 


19.45 


30.40 


13.54 


22.26 


17.30 


S~ 


32.00 


15.49 


16.07 


11.16 


17.97 


13.60 


n- 


19.79 


1.01 


8.03 


0.35 


1.17 


0.89 


n- 


19.79 


1.01 


5.53 


0.35 


1.17 


0.89 


7T + 


928.62 


911.72 


435.26 


1015.69 


1043.13 


811.42 




1877.62 


1843.35 


887.53 


2052.50 


2106.45 


1642.49 


7r~ 


948.99 


931.63 


452.26 


1036.51 


1063.32 


831.07 


K+ 


635.72 


686.80 


473.57 


630.45 


782.21 


613.97 


K- 


528.75 


570.90 


378.65 


534.03 


659.77 


503.06 


Kl 


589.21 


636.35 


435.36 


588.70 


728.55 


565.94 


Kl 


589.21 


636.35 


435.36 


588.70 


728.55 


565.94 


■n 


162.89 


125.6 


358.21 


86.23 


144.64 


111.06 



Table 6. The predicted hadronic observables in a central Pb + Pb collision at 5.52 TeV/nucleon-pair CM energy. 



16 is not farther from than \a. For the remaining 7 we 
get Table H 

There was really no measurement for the synthetic last 
bracketed ratio at 200 GeV RHIC, so that value is not 
analysed here. The last column is the error integral 

1 - erf{n) = ^/2/^ / e""'/^du (32) 

J n 

Were the distribution of A Gaussian, the last column 
would show the probability of such a random fluctuation; 
it is at least an estimation for it. The position is the posi- 
tion within the actual group in Table [71 

Individual measurements up to 2(7 are generally not 
regarded as too problematic. So 9 out of the 10 ratios 
outside Icr are not necessarily signals of problems, tests 
might be advisable, but will not be done here. However, 
the K~ /'K~ ratio at 200 GeV RHIC is surely problematic. 

Our guess is that we see consequences of enhanced pion 
yields, and it is probably enhanced by rehadronisation into 
resonances. As discussed already in Section 01 Point 2, 
the encounter of quarks may result in resonances instead 
of stables. Were the timescale of phase transition signif- 
icantly longer than the resonance lifetimes, the quarks 
could try again. However, both timescales are 10~^^ s, 
so the two processes do not clearly separate; some res- 
onances may decay in the hadronic phase, giving extra 
mesons, mainly tt's. 



We are not really manufacturing a full Model 12 in- 
cluding resonances; a lot of presently badly known new 
parameters should be used for that. But just for demon- 
stration consider the case when the final baryon yields 
are unchanged, but half of them come through the cor- 
responding member of the decuplet. (Of course, is an 
exception, being a stable and, in the same time, a mem- 
ber of the decuplet.) Then 3 model predictions of the 5 
200 GeV RHIC yields in Table [5] are unchanged, but 2 
gives new results. For estimation take the model errors 
equal to that of Model 7; then we get Table [HI That is our 
"Model 12"; The quotation marks are used because it not 
a fully elaborated model but rather an estimation. Both 
problematical ratios have got within Icr to 0. 

So surely there are stables which resulted through res- 
onances. Surely pion numbers will be then higher than in 
the present models. However, kaon numbers would be un- 
changed in first approximation since the decays mainly go 
as X* — >■ Xtt. So either we should handle the resonances 
properly before predictions are made; or as a first step 
we should ignore any ratio involving pions. As a second 
step, the TT surplus will be between and the sum of the 
numbers of 1/2-spin baryons -f antibaryons. 

The synthetic ratio K~ /p is still rather high from 
Model 7 compared to the ratios of 2 measured ratios. How- 
ever, the statistics of such derived quantities need serious 
cautions. 
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Energy /nucl. 


Ratio 


Experiment 


Reference 


Measurement 


Model 9 


158 GeV SPS 




WA97 


27 


0.14 ± 0.02 


U.14 ± U.13 




S'^/A 


WA97 


27 


0.26 ±0.05 


U.19±U.13 






WA97 


27 


0.19 ± 0.04 


U.27 ± U.13 






WA97 


27 


0.30 ± 0.09 


U.35 ± U.13 




A/ A 


WA97, NA49 


27 


0.145 ± 0.024 


U.38 ± U.13 




" /" 


WA97 


27 


0.27 ± 0.05 


U.49 ± U.13 




O/Q 


WA97 


27 


0.42 ± 0.12 


U.66 ± U.13 






AT \ AC\ 

IN A49 


27 


U.13 ± U.u3 


U.16 ± U.13 




/I In 




oo 


1 n^i + n 1 fi 


1 1 9 + n 1 

1 . _LZ ZIZ U.IO 






NA49 


28 


0.59 ± 0.05 


0.77 ± U.13 


zUU (jeV b-TO 


A/p 


IN Aoo 


Oft 
ZD 


U.oU d= U.ZO 


1 1 o _i_ n 1 /I 
i.iz zh U.i4 




/ A 
" //I 


WAoO 


1 o 


n 1 o -1- n ni 
U.ly =h U.Ul 


n 1 c 1 n 1 /I 
U. io ± U. 14 




/ A 
" //I 


WAoo 


1 O Oft 


n 01 -1- n no 








AT A o c: 
IN A^o 


26 


n on 1 n a 

U.oU ± U.4 


n OO 1 n 1 /I 

U.zo ± U.14 




/I / /I 


AT A Qc; 
IN Aoo 


Q Q 
OO 


n 1 Q _i_ n Oft 
U.iO ± U.UO 


U.oo ± U.14 












Model 7 


130 GeV RHIC 


viv 


STAR 


pi, [29] 


U.64 ± U.U5 


0.59 ± U.06 




Q/n 


STAR 


30 


1.UU±U.2 


1.00 ±0.06 




A/ A 


STAR 


28 , 29 


U.77±U.U5 


0.65 ±0.06 




j ^ 


STAR 


28 , 29_ 


U.81 ±U.U5 


0.78 ±0.06 






PHOBOS 


29_ 


U.95±U.U6 


1.U2±U.06 




R- /K+ 


STAR 


29 


U.9U±U.U5 


0.82 ±0.06 


200 GeV RHIC 


P/P 


PHOBOS, PHENIX 
STAR, BRAHMS 


.30,, 


U.84 ± U.U4 


0.60 ±0.16 




K- /K+ 


PHOBOS, PHENIX 
BRAHMS 


ED], IS] 


U.98 ± U.U4 


0.82 ±0.16 






PHOBOS 


3r 


l.U2±U.U2 


1.U2±U.16 




K- /%- 


STAR, BRAHMS 


30^ 


U.15±U.U2 


0.61 ±0.16 






BRAHMS 


30 


U.U8±U.U1 


U.16 ± U.16 




{K-/P 


STAR, BRAHMS 


30^ 


1.87±U.34 


3.81 ± U.16) 



Table 7. Experimental and model particle ratios for different collision energies. The errors of the penultime column are that 
of the ratio measurements. The average errors of the models can be obtained from Table [2] as: 0.13 for 158 GeV SPS, 0.14 for 
200 GeV SPS, 0.06 for 130 GeV RHIC, and 0.16 for 2UU GeV RHIC. 



Experiment 


Ratio 


A 


(7 


Deviation/cr 


Error integral 


158 GeV SPS 


A/A 


-0.235 


0.13 


1.85 


0.06 


158 GeV SPS 




-0.22 


0.14 


1.57 


0.12 


158 GeV SPS 


Q/Q 


-0.24 


0.18 


1.33 


0.18 


158 GeV SPS 


K- /K+ 


-0.18 


0.14 


1.28 


0.2U 


200 GeV SPS 


A/p 


0.32 


0.29 


I.IU 


U.27 


200 GeV SPS 


{Q + n)/{E- + E-) 


0.52 


0.42 


1.24 


U.22 


200 GeV SPS 


A/A 


-0.20 


0.15 


1.33 


U.18 


130 GeV RHIC 


A/A 


U.12 


0.08 


1.5U 


0.13 


200 GeV RHIC 


p/p 


U.24 


U.16 


1.5U 


0.13 


200 GeV RHIC 


K- /-K- 


-0.46 


U.16 


2.88 


0.004 



Table 8. 10 ratios in Table [7] whose a's are above 1 with one ratio whose ct's is above 2. The difference A of columns 5 & 6 is 
within Icr in 16 cases. 



Experiment 


Ratio 


A 


<j 


Deviation/(T 


2UU GeV RHIC 




-U.ll 


U.16 


U.69 


2UU GeV RHIC 


p/tv' 


U.U9 


U.16 


U.57 



Table 9. The so far problematic ratios at 200 GeV RHIC in "Model 12". 
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We emphasize that an almost infinite number of dif- 
ferent models can reproduce more or less the yields pub- 
lished so far. Our models are focused on rehadronisation, 
but they were not too interested about, say, collective phe- 
nomena in the QG or hadron phase. Such collective phe- 
nomena may cause quite remarkable effects, as e.g. [37] 
can predict very enhanced 77' yields, and then tt produc- 
tion, via the mass of 77' meson much lower than the free 
mass. Such effects are actually quite importans. However, 
note that 77', while it is technically a stable, is rather a 
short lived particle with 10~^^s lifetime and the situation 
is unfamiliar as in |16j shown many years ago that even 
existing particleis may or may not form an independent 
thermodynamic degree of freedom. In addition, lots of for- 
mally different approches may result in similar yields, as 
it was shown in [16^ . At the present status of knowledge 
we rather want to compare as many models as possible 
and we are not looking for the best one. 

At the end of this Section it is the proper place to 
note that there are very preliminary yields from the new 
range about Ecm = 2.76 TeV. E.g. [33] gives 2*4 ratios 
for nearly central collisions, without error bars. For an- 
tiparticles the ratios are almost exactly the same as for 
particles, and, according to what we mentioned above, it 
is better to leave out the pions. Then we remain with three 
ratios, /p, S~ /p and fi/p. Comparing these very ap- 
proximate ratios with predictions of Model 7 for 0.2 & 
5.52 TeVs we get Table [lOl 

While the first two ratios would be within error after 
the manner of Table 7, fi/p of Model 7 is too big, by cca. 
a factor of 5. However, Model 5 would give a ratio even 
bigger by almost 8. This again would suggest something 
between Models 5 & 7 as it was seen in Section [S] at RHIC 
energies. 



7 Conclusion 

At the end we should answer the question of QG phase. 
In some sense we have done that: we predict the yields 
of our Model 7, with S/N=65.7, Column 3 of Table U 
Now, Model 7 contains a phase transformation into QGP, 
then a definite amount of gluon fragmentation, and then 
a rehadronisation according to some rules. So if the yield 
ratios will be near to those calculable from our preferred 
Column 3, then indeed the detected ratios confirm the 
transient existence of the QGP. However, the old story of 
too much if+ in 1987 [51118115] is a warning that even a 
good agreement with the predicted ratios would not be a 
proof. A proof would be almost impossible. Observe that, 
e.g., a completely hadronic explanation of "too much" 
at 14.5 GeV fixed target could be "simply" an enhanced 
rate of the reaction 

p + p^p + K+ + A (33) 

above T=150 MeV. 

Namely, the existence of a QG phase is something 
with definite meaning only in thermodynamic context; 



and |21j gives the rigorous prescription to decide if a defi- 
nite thermodynamic phase is present. We cannot perform 
the needed operations in the later hadronic phase. In prin- 
ciple leptons may originate in QGP and so carry some 
direct information; but that is beyond the scope of the 
present paper. 

As told earlier, our preferred prediction is Column 3 
of Table [51 You can form then the actual ratios measured 
in the future. (Corrections for the slight d to u surplus at 
initial condition have been done in Table [51) The other 
columns are partly for comparison, partly for cases when 
"something unexpected happens" between TEVATRON 
and LHC energies. 

We have seen that at RHIC energies the performance 
of our Model 5 was only slightly worse that of Model 7; 
but combination of the two models did not help, so Model 
7 was chosen as the absolute best. Now, Model 5, with- 
out the selecting mechanism of hadronic final state data, 
results in more heavy hyperons than Model 7 (in Q the in- 
crease is almost 20- fold), and therefore, slightly less lighter 
ones. It is not quite impossible that with increasing energy 
the combinations feel less the final state masses; not be- 
cause of higher energy (entropy would not increase much) 
but because less time for new trials in hadron formation; 
but still we believe rather in Model 7. But data of [3S] sug- 
gest that the model might suppress fi too strongly, while 
Model 5 does not do it sufficiently. 

Model 11, Model 7 with enhanced gluon fragmenta- 
tion, would be observed, according to Column 5, as slight 
increases in lighter antihadrons and slight decreases in 
heavier ones. The differences are moderate, and it is ques- 
tionable if the observations will be able to distinguish 
Models 7 & 11. Yield ratios nearer to Model 11 than to 
Model 7 would not automatically suggest enhanced gluon 
fragmentation; rehadronisation into resonances and then 
decay in the hadronic matter leads to a tt surplus which is 
also true at enhanced gluon fragmentation. On the other 
hand, it seems that higher S /N (last column in Table [5]) 
would be distinguishable from resonance decay. Anyway, 
K~ /'n~ and p/tt", we think, again could signal resonance 
decays in hadronic phase if the ratios were much smaller 
than any columns of Table [5] "Model 12", as told, is not 
a full model but an estimation for the pion surplus from 
resonances. We simply assumed that during the hadro- 
nisation half of the emerging baryons were stables and 
the other halp resonances. The resonances decay in the 
hadronic phase. This "Model" has been explained at the 
end of Section [6] 

Column 6 shows the outputs for Model 7 at a higher 
S/N=75. The entropy would increase if something unex- 
pected happens with the cross sections between TEVA- 
TRON and LHC energies. For ratios the differences be- 
tween Columns 3 & 6 are slight, because all particle yields 
increase in Column 6. Maybe the difference would be ob- 
servable in the antiparticle/particle ratios; they are higher 
in Column 6, but the difference is not greater than 3 %. 

And finally, let us see the famous /'k^ ratio which 
provoked discussions back to 1987 about the formation 
of QGP EUmiS]; which lead to the ALCOR model in 
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Ratio 


0.2 TeV, Model 7 


2.76 TeV, prelimininary 
measurement 


5.52 TeV, Model 7 




2.8 ± 0.34 


3.2 


2.9 ± 0.49 


S-/P 


0.078 ± 0.013 


0.12 


0.083 ± 0.014 


n/p 


0.004 ± 0.0007 


0.02 


0.004 ± 0.0007 



Table 10. 

Comparison to non-pionic preliminary ratios for 2.76 TeV in [33. The errors of the model calculatios are based on Table [H] 





Year- 1987 


Model 5 


Model 7 


Model 9 


Model 11 


Model 7 


Note 


14.5 GeV fixed t. 


fragmentation 


favourite 


hadronic 


enh. fragmentation 


S/N=75 


K+/'K + 


0.24 


0.685 


0.753 


1.088 


0.621 


0.750 




^0 


0.557 


0.613 


0.837 


0.515 


0.620 



Table 11. Kj-K ratios in Pb+Pb reactions at 5.52 TeV. 



1999 [13] , and which was in early times in general expected 
a good QGP signal. For comparison we give K~ j'K~ in 
Table [m as well. 

Now, the morale is clear but somewhat disturbing. At 
"low energies" a Kj-n ratio much higher for kaons than 
for antikaons had been believed to be a quark signal. But 
for 5.52 TeV all models give higher j-K^ than K~ 
ratios; and the highest ratios are given by Model 9 com- 
pletely without QGP! (See also the K/ir discussion in 15|) 

One possible conclusion is that great care is needed 
when looking for signals. Another is that maybe at high 
energies the situation differs from those at lower ones. 
Both conclusions are indeed probable, even if not quite 
quantitative. However, there is something more, to which 
Table [TT] is a warning. 

Models 1-8 & 10-11 contain QGP, with quarks in equi- 
libria within (probably a good assumption, reaction times 
being guessed short in QGP); the rehadronisation assump- 
tions differ (and they may be not equilibrium ones). The 
completely hadronic Model 9 contains hadrons in equilib- 
ria. 

Now, for the QGP quark to antiquark surplus is neces- 
sary because of initial condition, except for s, where s/ s=l 
because all strange quarks were produced in QGP and 
strangeness is conserved. 

Being — {us) and — (us), the ratio K'^ / K~ > 
1 is nontrivial from quark composition. However, note that 
all models with QGP feel the initial conditions, so produ 
ce baryon to antibaryon surplus. Then /K~ > 1 is a 
necessary consequence of conservations. 

As for Model 9 without QGP it is an equilibrium model. 
If time were long enough, high K/t: ratios would be ob- 
tained. Now, in 1987 the time estimated seemed not enough 
for this, and we still believe this, because K~ /tt~ was cca. 
[IH]. Namely, consider initial N-l-N collisions. The sim- 
plest result is 3 particles, say N, a hyperon and a kaon. 
Antikaon instead of kaon is impossible for 3 particles be- 
cause of strangeness conservation. Therefore, reactions re- 
sulting in K~ would have low rates anyways; although the 
energy barrier is much lower in QGP for qqq — >■ qqqss then 
in HM for NN -> A^A'yl, cca. 300 MeV in the first case 
and almost 700 MeV in the second. So K~/K'^ << 1 in 



1987 in itself was not a signal either for or against liber- 
ated quarks, only against quarks in complete equilibrium; 
a "too high" K'^ number suggested liberated quarks, but 
because of K~ /K'^ << 1 is not a "quite proper" QGP. 

It seems that at RHIC energy, S/N=58 a proper QGP 
phase was formed. We cannot prove this being the QGP 
unreachable for us; but the hadronic yields prefer QGP. 
Then we have no reason to doubt in QGP at 5.52 TeV. 
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Appendix 



A List of the used rehadronisation models 

Here we use 12 models, that, however, are practically 8 
at LHC energies. They are as follows. Models from 1 to 8 
are rehadronisation model in 2"^ arrangement, as a binary 
number + 1 . The models were discussed in Ref. [S] & [13] ■ 
The conditions for the states of the 3 digits are as follows: 

1) Final state compressibilty is involved (1) or not (0). On 
this digit, however, dependence is negligible at 2.76-1-2.76 
TeV. 

2) Final state hadronic masses are taken into account 
(1) or not (0). This influences the yields and the best 
model cannot be selected a priori because the rehadro- 
nisation is not an equilibrium process. 

3) Gluon fragmentation takes place in QGP (1) or it does 
not (0). If it does, we take its parameters from Ref. j24) . 
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So Models 1-8 refer to the 3 digit states "000"-" 111" ac- 
cording to Table 1121 in summary. 



Model 


Gluon 


Hadronic 


Final state 


# 


fragmentation 


masses 


compressibility 


1 











2 








1 


3 





1 





4 





1 


1 


5 


1 








6 


1 





1 


7 


1 


1 





8 


1 


1 


1 



Table 12. A group of rehadronisation models including the 3 
options in columns 2-4. 



Model 9 is a purely hadronic phase one, and Model 10 
is a sequential fission one, see Refs. [5] & [T3]. 

Model 11 is a version of Model 7, but with gluon frag- 
mentation enhanced to the energetically possible maxi- 
mum, see Ref. dU] & [12. 

Model 12 is not an entirely elaborated model, but an 
educated guess, when everything starts as in Model 7, but 
half of the resulting baryons enter the hadronic phase as 
resosances and the decay there, giving surplus mesons, 
mainly tt's. 



B Transformation to partially specific 
variables 

Eq. (1121114^ guarantee that we can calculate functions from 
s. Now, let us use s. 

As told Uneq. ([S]) must hold, for any timelike field. 
In s-r the derivatives n^ru'' and e^rU^ can be substituted 
from Eqs. (|3I4I) and then we get 

S,t Dt + (J - ns,n -es,e -ks,e )u'';r +t''{z,r ~S,e /3,r ) 
> 

D = u'ds (B.l) 

identically i.e. for any u; and several terms apppear with 
quite various u-dependences. Collecting the terms with 
the same u-dependences, Uneq. ([5]) is guaranteed only if 
the terms of different u-dependences satisfy equations in- 
dependently. That is set in Eq. [15]) . 



C Evaluation of A of ultrarelativistic collisions 

This Appendix is necessary to avoid inconsistencies in the 
main text. Namely, in this paper t stands for 3 quite dif- 
ferent quantities: time, the specific extra extensive and 
one invariant familiar in 2-particle collisions. Similarly, s 
means both entropy density and another invariant, while 
u is both the velocity vector and the third invariant. So 



we have separated the part most confusing and here we 
try to be careful and explicit in the formulation. 

Consider an elastic two-body collision where the in- 
coming two particles have the original momenta of the 
beams. By appropriate choice of the coordinate system 
and by momentum conservation we always can write: 

m1/2 = {(1 + z;2)1/2^±„^0,0} 

«3/4 - {(1 + V'y/^±{V^ - WY^^±W, 0} (C.l) 

Now, the collision is characterised by a triad (s,t,u) which 
here for clarity will be written in bold. They are defined 
as 

t = (Pi - Psf 

u={p,-p,f (C.2) 

but here we will not use u. Using Eqs. (jC.llC.2l) we get 

s = 4m2(l -f (C.3) 
t ^ -2m^v{v - {v^ ~ w^)^/'^) (C.4) 

Expressing it with the decrease of v, (u^ — w^) = v — A, it 
is simply 

A = t/2m^c^v (C.5) 

The average loss of the longitudinal momentum per colli- 
sion is then 

< A>= {l/2Tn^c^v){ J tcr(t)dt/ J cr(t)dt) (C.6) 

where v can be substituted by s via Eq. (|C.6p . 

Since A is the source term of q, the density of the 
momentum-like extra extensive, A must be the product 
of 3 terms: the collisions/time calculated from the total 
cross section, the average momentum loss in one collision 
(this is < A >), and the actual density. The last is cca. 
2no. As for the s dependence we used the almost Iri^s 
curve of and [55]. 

D The A-E sector 

Among the stables we find 4 hyperons with (qqs) quark 
composition: (uus) is , (dds) is S~ , but (uds) may 
both be 17° and A. The difference is in the relative spin 
positions. Measurements indicate that (qq) is parallel in 
S and antiparallel in A; and the antiparallel configuration 
is absent for {uu) and (dd) (For S is 1(5^), for A 

is 0(^^) ). No doubt, we see a consequence of the Pauli 
principle. 

Of course, in the spin-3/2 decuplet antiparallel spin 
position is per definitionem impossible, so there only a 
triplet S* appears with — 1(|^). The interesting fact 
is the 77 MeV mass difference between A and A is 
preferred. Obviously this comes from the spin positions. 
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but our rchadronisation riiodcils nciglcctcd the energy dif- 
ference within the sectors. But then the yields amongst 
the Y's are debatable from purely the quark numbers, 
while the total Y yield is unique, and the differences do 
not influence the other baryons. We think that the final 
answer would need the handling of resonances in rehadro- 
nisation, and in this moment that would delay the predic- 
tions. So instead we elaborated 2 simple models and chose 
that which kept charge conservation better. Here we dis- 
cuss only the charge-symmetric distinction of A and S; 
afterwards the inclusion of u/d ^ 1 is as for the other 
particles. 

Proposal 1 was S=2*Y/3, A=Y/3. This assumes that 
every initial spin position remains if a stable is possible, 
and fades away if not. So is Y/6 and = A =Y/3. 

Proposal 2 is = yl =Y/2, the parallel and antiparal- 
lel positions of the non-strange quarks are equally proba- 
ble. Then the charged pre-yl's first form E*'s, which then 
decay into A and tt's. 

While Proposal 1 seemed more logical, it generally re- 
sulted in Z « 174, while Proposal 2 in Z 172. Exact 
charge conservation would be Z=164. Probably, A > 
because the A mass is less; but this would not help the Z 
conservation, and here for now we stop. 
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